%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Calcul de la dérivée seconde en un noeud autre que point limite
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function [deriv0, derivl] = Derive_seconde_partielle(uhloc,l)

deriv0 = -6/(l^2)*uhloc(1)-4/l*uhloc(2)+6/(l^2)*uhloc(3)-2/l*uhloc(4);
derivl = -6/(l^2)*uhloc(1)-4/l*uhloc(2)+6/(l^2)*uhloc(3)-2/l*uhloc(4)+...
         +12/(l^2)*uhloc(1)+6/l*uhloc(2)-12/(l^2)*uhloc(3)+6/l*uhloc(4);

end